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We show that the main results of the analysis of the friction factor for turbulent 
pipe flow reported in G. Gioia and P. Chakraborty (GC), Phys. Rev. Lett. 96, 
044502 (2006) can be recovered by assuming the Heisenberg closure hypothesis for 
the turbulent spectrum. This highlights the structural features of the turbulent 
spectrum underlying GC's analysis. 
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INTRODUCTION 



The accurate prediction of the friction factor for the turbulent flow of an incompressible 
fluid through a straight pipe of constant circular section is a matter of huge importance, 
both practical jmd fundamental Q, [2], 3}. Not surprisingly, it has been the subject of careful 



measurement 



in concrete situations 



4j] . There is also an array of empirical formulae to predict the friction factor 

nn 

[5|, |6J , culminating with the sheer rendition of Nikuradse's experimental 

7]. However, the theoretical link 



results in analytical form provided by Yang and Joseph 
between the phenomenological formulae and the experimental results is weak. In particular, 
some of the most used empirical formulae, such as Colebrook's, erase most of the structure 
seen in experiment. 

In these circumstances, the derivation of several of the key features of the dependence of 
the friction factor with respect to Reynolds number for a given pipe roughness from a concrete 
theoretical model in 8fl (henceforth called GC) is undoubtedly an important step forward. 

nnnnn 

For some background and further developments on the GC model see [10, Hll, H2, 1131, 1141 ]. 
The analysis in GC departs in important aspects from the view of the problem laid down 
by pioneers such as Prandtl and von Karman [1] and the classic textbook formulation by 
Landau |9|]. Therefore it is important to understand what are the fundamental elements 
underlying the success of the GC model. 

Several key features of the Reynolds number dependence of the friction factor, namely the 
bellies and the Strickler's regime (see GC), can be regained if the GC analysis is combined 



with Heisenberg's closure hypothesis 
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171 ] for the turbulent spectrum. The Heisenber 



theory is not generally regarded as a realistic depiction of fully developed turbulence |18l.ll9l]. 
Therefore, the fact that the GC analysis works even if the at best qualitatively correct 
Heisenberg theory is used instead of a (not yet known) exact turbulent spectrum gives us a 
new perspective on the inner working of the GC model. 

This paper is organized as follows. The next section is of a review nature and presents 
the basic facts and definitions concerning the friction factor and the GC treatment thereof. 
Section 3 presents the Heisenberg closure hypothesis and derives the GC friction factor for 
this form of the spectrum. We conclude with a final appraisal of the GC friction factor 
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formula. 



II. PRELIMINARIES 



A. Definitions 



The goal of this Section is just to put together a basic theoretical description of the 
Reynolds number dependence of the friction factor. Therefore we shall consider only the 
simplest case of a single phase incompressible fluid moving within a horizontal pipe. The 
section of the pipe is circular and the radius is R. We assume the flow is well developed, 
meaning that there is a well defined macroscopic velocity V at every point (we shall use 
boldface for vector quantities, except for the Reynolds number Re to be defined below). V 
depends only on the radial coordinate r, points in the axial direction x, and vanishes at 
the boundary: V = v (r) x and v (R) =0. V is automatically divergenceless and obeys the 
momentum balance equations, which, written in cylindrical coordinates, are 

(1) 

ox r or 
Or r or 

Where p is the pressure and r is the stress tensor. Since the RHS of eq. [2] depends only on 
r, we the pressure drop dp/dx must be r- independent. The eq. [I] can be integrated (with a 
boundary condition imposed by regularity) 

rx r dp _ f r \ 

where r is the stress at the wall. This means that to know the stress on the wall it is enough 
to find the stress anywhere, since it obeys a simple scaling law. 

Another important quantity is the average flow Q. Together with the mass density p and 
the cross surface A = ttR 2 it defines the mean velocity V according to Q — pAV. The fluid is 
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characterized by a molecular (dynamic) shear viscosity p, (kinematic) v — pj p. With these 
quantities we may construct the most important dimensionless number, namely Reynolds' 

2W 

Re = (4) 

v 

These scales allow us to construct an energy density scale e = p V 2 . The Darcy-Weisbach 
formula introduces the friction factor / from the ratio of tq to e 

to = {p V 2 (5) 
Our main goal is to find a relationship between / and Re. 



B. Main flow regimes 

In this subsection we shall describe the main flow regimes and the corresponding empirical 
formulae. We shall use as reference Nikuradse's experimental results. For practical reasons, 
we do not mean the actual results, but rather Yang and Joseph's analytical rendering thereof 

For a rough pipe the friction factor does not decrease indefinitely with increasing Reynolds 
number, but rather converges to a finite value /oo. This allows us to define a parameter 
5oo = eR from the condition that 

(6) 

We shall use this parameter to identify the several series of data from the Nikuradse exper- 
iment . 

Fig. [1] gives an overall impression of the data. Each curve shows a rich structure. For 
example, let us consider the curve corresponding to e _1 = 126 (Fig. [2]). The log-log plot 
is essentially linear to the left of A. This corresponds to the laminar flow; at A, we have 
Re ~ 2, 000. Then there is a maximum (the so-called hump) at B, corresponding to Re ~ 
4,000. The log-log plot is again essentially linear up to the belly at C (Re ~ 40,000). To 
the right of C, the curve approaches its asymptotic value from below. The approach is very 
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FIG. 1: Analytic reconstruction of Nikuradse's data, as given in [7J, extrapolated to the range up 
to Re = 10 6 , actually covered in the Princeton Super Pipe experiment. The six lines correspond, 
from the bottom up, to e _1 = 507, 252, 126, 60, 30.6 and 15. Both scales are logarithmic. 



fast; beyond the point D (Re ~ 200, 000) the friction factor is constant for all practical 
purposes [7] 



0.05 



0.02 
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FIG. 2: A close up of the curve corresponding to e 1 = 126. Observe the sharp minimum at A, the 
hump at B and the belly at C. To the right of D the friction factor is constant. 



We shall now review two basic regimes in this complex behavior, namely, the laminar 
regime, to the left of A, and the Blasius regime, from the hump B to the belly C. 
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1. The laminar regime 

If the flow is laminar, the bulk velocity profile can be solved exactly and the friction factor 
results 




FIG. 3: A comparison of the friction factor for laminar flow eq. [7] (thick line) against Nikuradse's 
data as given in fig. [1] 

In fig. [3]we superimpose the plot of the friction factor for laminar flow eq. [7]to Nikuradse's 
data as given in fig. [TJ We can see that the agreement is outstanding up to Reynolds numbers 
of a few thousands. 



2. The Blasius regime 



One of the oldest and most accurate empirical formulae for the friction factor is Blasius' 

f - (R) 

J Blasius — V ' 

We plot this expression superimposed to Nikuradse's data in fig. HI 
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FIG. 4: A comparison of Blasius' friction factor eq. [8] (thick line) against Nikuradse's data as given 
in fig. [TJ 



C. The GC model 

In this subsection we shall analyze the proposal by Gioia and Chakraborty (GC) 8]. Our 
aim is to contrast it with the experimental situation as described above. 

The basic framework of the GC model is that for high enough Reynolds number the 
pipe is filled with a well developed turbulent flow which can be accurately described by 
Kolmogorov's K41 theory. It is therefore characterized by large velocity and length scales 
Uqc and L GG . On scales I smaller than L GG , the turbulent speed follows Kolmogorov 
scaling uqc (0 = Ugc (I I 'Lgc) 1 ^ (except for small corrections to be discussed below). GC 
identify L GG — R as the pipe radius. Concerning the macroscopic flow velocity, GC assume 
Ugc — K vV i where V is the mean flow velocity, and K u IS cl He- independent constant 
(eventually GC choose k u = 0.036). 

The detailed mechanisms of momentum transfer between the flow and the wall are rather 



complex 20]. GC assume that the transfer is mainly effected by eddies of size 5qc- These 
eddies carry a momentum pUcc along the wall. The transverse velocity u° GC is the turbulent 
speed associated with the scale Sqc 



u gc ( S gc) = U GC (5 GC /L GC ) 1/3 



(9) 



The stress at the wall is r = k t pUgc u gc with some universal constant k t . Since To is linear 
on u GC , in the Blasius regime 5 GC oc Re _3/4 . This is the same scaling as Kolmogorov's 
viscous length scale rj. 77 = i?/Re -3 ^ 4 is the scale at which viscous dissipation matches 
turbulent energy transport (GC actually interpose another dimensionless constant b = 11.4 
in the definition of 77, we prefer to absorb it into the many more constants to come). Therefore 
it becomes natural to assume 5gc — ar ]i where a is yet another universal dimensionless 
constant. The friction factor reads 

f = ^J^ F ^ ( 10 ) 

where 

At the level of approximation we have stayed so far, F = 1. If we choose 8K T K^a^ 3 = 
0.3164, then the model is built to reproduce Blasius' law for smooth pipes. This can be 
criticized on the grounds that it does not account for the deviations from the Blasius' Law 
at very high Reynolds number predicted by the Colebrook equation and apparently verified 



by the Princeton Superpipe [4] . However, these criticisms can be swept away by noting that 
no physical pipe can be absolutely smooth. So the real issue is how to introduce roughness 
into the model. 

The way to modify the asymptotic behavior of the friction factor to account for pipe 

roughness is that the scale Sqc does not decrease indefinitely with growing Reynolds number, 

1/3 

but rather stabilizes at a value ecc'R- We therefore find /oo oc 6q C . 

If we identify ecc with the parameter e in eq. El then this result is the so-called Strickler 
Law. This approximate scaling was known long before GCs work, and the fact that it is so 
effortlessly obtained is quite remarkable. 
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1. The belly and the Blasius regime 

For finite Reynolds number, we may imagine that Sqc — ^gc-R + a V- This is, the width of 
the dominant eddies at finite Reynolds number is just the sum of the widths defined by the 
pipe roughness and by the viscous scale. However, if the stress at the wall is defined by the 
turbulent velocity at Sac, and this in turn decays to eccR with increasing Reynolds number, 
then Kolmogorov scaling implies that the friction factor reaches its asymptotic value from 
above, as in the Colebrook equation [||. To recover the belly in the Nikuradse's data (cfr. 
figs. [H and [2]) we need that the velocity at an essentially fixed scale e^c-R be an increasing 
function of Re, at least for large enough Reynolds number. 

The basic idea is that the mean square velocity at some scale 5 is the sum of contributions 
from eddies at all scales k~ l < 5 

roo 

u 2 [S] = / alkS[k] (12) 



/ is proportional to (u 2 [S]) 1 ^ 2 . If the only dependence on Reynolds number were through S, 
then we would obtain 

it _ s «n <« <0 (13) 



dRe 2u 2 [5] 5 2 dRe 
Therefore the friction factor would be a monotonic function of Reynolds number. To get 
the feature of the belly in the friction factor we have to allow for a direct dependence of the 
spectrum on Reynolds number, leading to the proper result 

if _ f (r ft «H + ™ (14) 



dRe 2u 2 [5] [J 6-i dRe 5 2 dRe 
The derivative does not have a definite sign. 

As a matter of fact, this sort of behavior is a prediction of K41 theory for finite Reynolds 
number. For large but finite Reynolds number, we have 



(15) 
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where E is a non increasing function such that E [0] = 1 and E [oo] = 0. Thus in the limit of 
infinite Reynolds number i^nOwe recover the one-third scaling law eq. [91 We have added 
the (yet another) dimensionless constant (3 for latter convenience. Neither f3 nor the form of 
E more generally are prescribed by K41 theory. This is precisely the point where adoption 
of the Heisenberg closure hypothesis (or any other, see 22]) makes a difference. 

We observe that inclusion of the E \f3krj\ factor enforces a decay faster than —5/3 in the 
dissipative range. In any case, we now have 



dS[k] 



3 3U GC 



1 (3 V 



-E' \pkrj\) > 



(16) 



dRe 4 2i?2/3 fc 2 /3Re 
So the derivative in IT41 may be positive in the appropriate range. 

We can derive a more explicit result. Using eq. IT51 into 1121 and calling x = fikr] we get 
after an integration by parts, 
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S) 



2/3 



dx 

/3rj rr>2i 1 3 

T ^ 



{-Ef [x\) 



(17) 



If we use eq. 1171 to evaluate u at the scale 5gc = ^GcR+^r], we clearly see the two asymptotic 
regimes. For large Reynolds number, 77 1 — >• 0, 5gc ^ £gcR and u = Uqc^gc or ' retaining 
first order corrections 



u = Ugc^gc \ 1 



B 



provided E [(3t]/5] = l + o [r] 2 ^ 3 j, where the constant B depends on the shape of the function 
E as well as on a, ecc and j3 



2/3 



dx 



2\e GC ) Jo x^ 
In the limit where ar] ^> eccR w e get instead 



-E'[x\) 



(19) 



u 2 [av] = U GC (^) 2/3 \e 



' j3\ 2/3 r°° dx 



-E'[x\) 



(20) 



a I J ( - x 2 / 3 

/ a 

and so we recover Blasius' Law with a new constant (unless /3/a <C 1). Of course this new 
factor can be easily absorbed in any of the several dimensionless constants at our disposal. 
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This qualitative success is not easily transformed into a quantitative fit, however. Gener- 
ally speaking, to get the fast approach to the asymptotic value characteristic of Nikuradse's 
data, very small values of (3 are preferred. But then the Blasius regime appears in a range 
of Reynolds numbers much below the experimentally observed. In any case, we have not yet 
accounted for the hump. We shall defer further quantitative analysis until we incorporate 
the hump into the GC model. 



As we have seen, the behavior of the friction factor for Reynolds numbers over a few 
thousands, according to GC, is the result of the competition of two opposite processes. On 
one hand, higher Reynolds numbers mean higher lower limits in the integral eq. [121 thus 
bringing the friction factor down. On the other, the integrand in eq. O as defined in eq. 
[T5l increases pointwise with Reynolds number, thus bringing the friction factor up. The first 
process dominates in the Blasius regime, the second in the climb up to the asymptotic value. 

To set a lower limit to the Blasius region, therefore, the simplest is to cut off the integral 
in [TJl so that it becomes insensitive to the lower limit if this is low enough. The first process 
then becomes moot, while the second is still operative, and we get a friction factor which 
grows with Reynolds number. 

We may mention that this second modification of the spectrum is totally outside K41 
theory. Also that in a certain way it works too well, since the laminar regime is obliterated. 

In summary, GC propose the form of the spectrum (cfr. eq. fT5i) 



2. The hump 



£[k] 



2R2/3 k 5/3 



E \j3krj\ D [kR] 



(21) 



where 



E[x) 



= e 



—x 



(22) 



x 



17/3 



D[x] 



[x 2 + 7] 



17/6 



(23) 
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and 7 is a dimensionless constant. Introducing a dimensionless integration variable x = (3krj 
the friction factor may be reduced to the form 

C ( f°° x 4 e~ x \ 1/2 

where 

h [ Re l = 7 \ 3/4 ( 25 ) 



5 [^] = (26) 



7^ 2 
Re 3 

In fig. [5]we show a typical plot of eq. [24j corresponding to the data for e = 1/126. We have 
extrapolated the experimental data up to Re = 10 10 to better appreciate the convergence to 
the asymptotic value. For this plot, we have chosen f3/a = 0.5 and 7/3 2 = 10 4 . The values of 
the constant C and of eqc were chosen to enforce the Blasius law at the value of Re = 4000 
and the proper asymptotic limit. 

This plot represents a near optimal situation; for example, higher values of 7 erase the 
features of the curve, while for lower values the hump disappears. Although there is a 
resemblance to the experimental data, it is not truly quantitatively accurate; in particular, 
the theoretical prediction is much smoother than the experimental plot. 

III. A GC MODEL BASED ON HEISENBERG CLOSURE 

In this section we shall use the Heisenberg closure to derive the form of the spectrum in 
eq. [121 thereby obtaining an expression for the friction factor with higher predictive power 
than the original GC proposal. We consider this as a toy model to test whether a turbulent 
spectrum based on a definite hypothesis regarding the underlying dynamics can make a 
difference in the accuracy of a model built along the general lines of the GC approach 23]. 
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FIG. 5: Log-Log plot of the friction factor from eq. [241 compared to the data for e = 1/126. 
The thick line corresponds to the GC model, the thin line to the data. We have extrapolated the 
experimental data down to Re = 10 2 and up to Re = 10 10 to better appreciate the convergence to 
the asymptotic values. For this plot, we have chosen f3/a = 0.5 and 7/3 2 = 10 4 . The values of the 
constant C and of ecc were chosen to enforce the Blasius law at the value of Re = 4000 and the 
proper asymptotic limit. 



The basic idea of Heisenberg closure is in fully developed turbulent flow the sum of the 
energy dissipated by viscosity and the energy transported to smaller eddies remains the same 
in all scales. At a given scale k , the energy dissipated through viscosity is 



A. The Heisenberg closure in the presence of boundaries 




(27) 



where 





Heisenberg assumes that the energy transported by turbulence may be written as 




(29) 
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where, on dimensional grounds 



DC 



Vturb [k ] = A ^ dk^^- (30) 
where A is a dimensionless constant. The basic assumption is thus 

Qu [h ] + Qturb [k ] = Qtot (31) 

independent of k . If £ [k] oc k a when k \— > 0, then (Vw 2 ) [k ] scales as k$ +3 while v turb 
scales as k^ 2 1//2 , so we must have a = —5/3. This means that the kinetic energy increases 



indefinitely with eddy size, which is incompati 
adopt the modification suggested by Parker 
and M by 



ale with the presence of the pipe. We therefore 
24j . namely, we replace (Vu 2 ) [ko] in eqs. E7] 



y u 2 ) R [k ] = p dk (k 2 + K 2 ) 8 [k\ (32) 

where K « The new term accounts for the increase in dissipation due to eddy defor- 
mation. Our power counting argument now gives a = —1/3, so that the total kinetic energy 
in the flow is finite. 

The starting point is than the balance equation 



{u + u turb [k] } 2 ( Vu 2 ) r [k] = Qtot (33) 



A derivative of eq. [331 gives 



but also 



#^(V^ 2 ) ^~ 2 VS = (34) 

{Vu 2 ) R dk\ ' R \ ' R k z l 2 v 1 



£ = —±—4r (v« 2 ) (35) 

k 2 + K 2 dk\ IR v ' 



so, if £ = k 2 
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we read the boundary condition off eq. ED 



A -3 



R 



Q 3 



+ 



6A 



tot 



K 2 Q 



2 

tot 



k 2 K 2 ln 



K 2 



(37) 



and then, from eq. [35], we get the spectrum from a simple derivative 



£ 



1 



AA [ 8z/ 3 



k 3 (k 2 + K 2 Y Qtt 
We adopt units where K = R^ 1 
Introducing 



Qtot 



+ 



6A 



K 2 Q 2 tot 



1 J-l 

k 2 K 2 n 



4/3 



1. In these units, Q 



lot 



k 2 

V 3 and v 



3/4 



(38) 
2y/Re. 

(39) 



then the spectrum can be rewritten as 



£ = constant — - 



3 1 ^ + 2(F-" 1 



1 + F 



-4/3 



3fc 3 (P + l) 2 K s 4 ' " ' " ' 1 ' IU) 

To obtain the friction factor, we need to integrate the spectrum from a lower scale Kq up 

to infinity, where 



(41) 



{^HeiK s + ClHei) 

After identifying K s = 1/rj, this lower limit is the same as in the GC approach. 

Observe that the spectrum l4*0l behaves as A; -5 / 3 in the inertial range. It displays a faster 



.a 



decay Re 4 /c -7 in the dissipative range. It therefore satisfies the criteria discussed in ref 
and in the previous Section for reproducing the Blasius and Strickler scaling, as well as the 
belly feature in the friction factor plot. 

To compare the performance of the model based on Heisenberg closure with the results 
from the original GC proposal, we shall seek an optimal fit to the Nikuradse data for e = 
1/126. In all, we have four parameters Re c , ciHei, ^Hei and an overall normalization at our 
disposal. To reduce parameter space, we shall assume esei — 1/126 as well. The overall 
constant is determined by asking for a good fit for very large Reynolds number. Changes in 
Re c induce rigid horizontal shifts in the plot, so a is the parameter which controls the shape 
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of the curve. The best fit is obtained for a = 1.25, which is physically acceptable. We show 
the result in fig. [6] 

f 

0.05 
0.02 

10 2 10 3 10 4 10 5 10 6 

FIG. 6: Log-Log plot of the friction factor as computed from the form[3U]of the spectrum, compared 
to the data for e = 1/126. The Heisenberg model prediction is the thick line, the experimental 
values are the thin line. We have extrapolated the experimental data down to Re = 10 2 and up 
to Re = 10 10 to better appreciate the convergence to the asymptotic value. For this plot, we have 
chosen a = 1.25, eHei = 1/126 and Re c = 1000. The overall constant is determined by fitting the 
theoretical value for / at infinitely large Reynolds number to the experimental asymptotic value. 

We see that the Heisenberg closure leads to an expression which is as successful as the 
GC model in describing the Blasius regime and the approach to the asymptotic limit. This 
is clearly displayed in fig. [3, which is simply the superposition of figs. [5] and [61 

IV. FINAL REMARKS 

The GC proposal is striking in that it offers a simple explanation for the ups and downs 
of the plot of the friction factor against Reynolds number. 

In this note we have combined the main GC postulates with a spectrum derived from 
Heisenberg closure. The result fits the Nikuradse data as well as the original GC analysis. 
This is good news for the GC model in that it underlines the fact that the model is built on 
generic features, rather than detailed dynamical characteristics, of the flow. However, it is 
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FIG. 7: The plots of figs. [5] and [6] combined. The experimental data corresponds to the short 
dashes, the GC model is the long dashes, and the Heisenberg closure is the full line. There is 
essentially no difference between both theoretical models up to the highest Reynolds numbers. 

also problematic, because Heisenberg closure is not generally regarded as realistic, specially 
in the dissipative range. 

Both in the original and the Heisenberg closure model, moreover, it is clear that the 
quality of the final fit depends on the careful tuning of the many available parameters. 

Our conclusion is that to make progress in understanding the friction factor of turbulent 
pipe flow along the direction pioneered by Gioia and Chakraborty we do not need to worry 
about a more accurate spectrum shape, but rather to provide a solid foundation for the basic 
assumptions of the GC analysis. 



We wish to thank Jose Lanziani, Salvador Gil and Pablo Mininni for discussions and 
pointing ref. j^] to us. 

This work is supported by University of Buenos Aires, CONICET and ANPCyT. 



Acknowledgments 



[1] H. Schlichting, Boundary Layer Theory (McGraw-Hill, New York, 1950). 



18 



[2] A. S. Monin and A. M. Yaglom, Statistical Fluid Mechanics (MIT Press, Boston, 1971). 
[3] S. B. Pope, Turbulent Flows (Cambridge UP, London, 2000). 

[4] B. J. McKeon, M. V. Zagarola and A. J. Smits, A new friction factor relationship for fully 

developed pipe flow, J. Fluid Mec. 538, 429 (2005) 
[5] G. O. Cordero, Improvements in Estimating Pipeline Friction, GIMOR (SPE International- 
Argentine Petroleum Section-Grupo de inters en modelado y operacin de redes y ductos) 7th 

Annual Meeting, 5 November 2008, Buenos Aires. 
[6] N. Afzal and A. Seena, Alternate scales for turbulent flow in transitional rough pipes: universal 

log laws, Transactions of the ASME 129, 80 (2007). 
[7] B. H. Yang and D. Joseph, Virtual Nikuradse, J. Fluid Mech. (submitted) 

http://www.aem.umn.edu/People/faculty/joseph/PL-correlators/docs-ln/Rl-JFM-submissions-virtual-Nikui 
[8] G. Gioia and P. Chakraborty, Turbulent friction in rough pipes and the energy spectrum of 

the phenomenological theory, Phys. Rev. Lett. 96, 044502 (2006). 
[9] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Pergamon Press, Oxford, 1987) 
[10] G. Gioia and F. A. Bombardelli, Scaling and similarity in rough channel flows, Phys. Rev. 

Lett. 88, 014501 (2001). 

[11] G. Gioia, P. Chakraborty and F. A. Bombardelli, Rough-pipe flows and the existence of fully 

developed turbulence, Phys. Fluids 18, 038107 (2006). 
[12] N. Goldenfeld, Roughness-induced critical phenomena in a turbulent flow, Phys. Rev. Lett 96, 

044503 (2006). 

[13] M. Mehrafarin and N. Pourtolani, Imtermittency and rough-pipe turbulence, Phys. Rev. E77, 
055304(R) (2008) 

[14] N. Guttemberg and N. Goldenfeld, The friction factor of two-dimensional rough-pipe turbulent 
flows, larXiv:0808.145Tl 

[15] W. Heisenberg, On the theory of statistical and isotropic turbulence, Proc. Roy. Soc. A195, 
402 (1948). 

[16] S. Chandrasekhar, On Heisenberg's elementary theory of turbulence, Proc. Roy. Soc. 200, 20 
(1949). 



[17] E. Calzetta and M. Grana, Reheating and turbulence, Phys. Rev. D65, 063522 (2002). 
[18] G. K. Batchelor, The theory of homogeneous turbulence (Cambridge University Press, Cam- 
bridge (England), 1953) 
[19] J. O. Hinze, Turbulence (McGraw-Hill, New York, 1975). 

[20] A. A. Townsend, The structure of turbulent shear flow, (Cambridge University Press, London, 
1976). 

[21] A. Yaglom, The century of turbulence theory, in New trends in turbulence theory, Les Houches 

course LXXIV, edited by M. Lesieur, A. Yaglom and F. David (Springer Verlag, Berlin, 2001) 
[22] V. L'vov, I. Procaccia and O. Rudenko, Universal model of finite Reynolds number turbulent 

flow in channels and pipes, Phys. Rev. Lett 100, 054504 (2008). 
[23] C-M Tchen, Transport processes as foundations of the Heisenberg and Obukhoff theories of 

turbulence, Phys. Rev. 93, 4 (1954). 
[24] E. Parker, Extension of Heisenber's model of turbulence to critical Reynolds numbers, Phys. 

Rev. 90, 221 (1953) 



